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Abstract. We study the one-dimensional Holstein model of spinless fermions interacting with dispersion- 
less phonons by using a recently developed projector-based renormalization method (PRM). At half-filling 
the system shows a metal-insulator transition to a Peierls distorted state at a critical electron-phonon 
coupling where both phases are described within the same theoretical framework. The transition is accom- 
panied by a phonon softening at the Brillouin zone boundary and a gap in the electronic spectrum. For 
different filling, the phonon softening appears away from the Brillouin zone boundary and thus reflects a 
different type of broken symmetry state. 
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1 Introduction 

Systems with strong electron-phonon (EP) interactions 
have received considerable attention in the last few years, 
motivated by new findings which suggest a crucial role of 
the EP coupling in materials with strong electronic corre- 
lations. Examples are high-temperature superconductors 
PP, colossal magnetoresistive manganites 0, or metallic 
alkaline-doped Cgo-based compounds [3]. Furthermore, in 
many quasi one-dimensional materials, such as MX chains, 
conjugated polymers or organic transfer complexes [3], 
the kinetic energy of the electrons strongly competes with 
the EP interaction which tends to establish, e.g., charge- 
density wave structures. 

The so-called spinless Holstein model, 

n = - l E( c 1 c j + h - c -) + ^ E ^ w 

+.?E {b\+h) ni , 

i 

is perhaps the simplest realization of a strongly coupled 
EP system. It describes the local interaction g at a given 
lattice site i between the density m = c\ci of electrons 
and dispersion- less phonons with frequency ljq. Here, the 
c\ (b\) are fermionic (bosonic) creation operators of elec- 
trons (phonons), and denotes the summation over all 
neighboring lattice sites i and j. 

In the past, a large number of different analytical and 
numerical methods have been applied to the Holstein mo- 
del . Strong-coupling expansions |S| , variational [B| and 



renormalization group 0E| approaches, as well as Monte 
Carlo [S]0|] simulations were used to investigate mainly 
ground-state properties. More recently, exact diagonaliza- 
tion [10] (ED) and density matrix renormalization group 
(DMRG) jl lU12U13j techniques, and dynamical mean-field 
theory (DMFT) in conjunction with a numerical renor- 
malization group approach |14j were applied. However, 
most of these approaches are restricted in their applica- 
tions, e.g., ED techniques can only handle very small sys- 
tem sizes, DMRG methods require one-dimensional sys- 
tems, and the DMFT exploits the limit of infinite spatial 
dimensions. Furthermore, the phononic part of the Hilbert 
space is infinite even for finite systems so that all numer- 
ical approaches require truncation schemes to limit the 
number of bosonic degrees of freedom, or a numerically 
expensive systematic reduction of the Hilbert space in the 
spirit of the DMRG method JS] has to be employed. For 
these reasons, there is still a clear need of reliable theo- 
retical methods to tackle strongly coupled EP systems in 
terms of minimal theoretical models. The Holstein model 
of spinless fermions shows at half-filling a quantum- 
phase transition from a metallic to an insulating state 
where both the one-dimensional case and the limit of in- 
finite dimensions have been studied [olll 11114) . 

Alternative analytical approaches to interacting many- 
particle systems are offered by methods based on func- 
tional renormalization like the flow equation method |16j , 
the similarity transformation |17| . or the projector-based 
renormalization method (PRM) JSj. So far, these meth- 
ods have been successfully applied to electron-phonon sys- 
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tems with the aim to study the effective electron-electron 
interaction or superconductivity |2()| . However, the 
quantum-phase transition in the Holstein model has not 
yet been studied by this kind of approach. 

Recently, we applied the projector-based renormaliza- 
tion method (PRM) ^2] to the spinless Holstein model 
l|T]) at half- filling [21| . This analytical approach does not 
suffer from the truncation of the phononic Hilbert space 
so that all bosonic degrees of freedom are taken into ac- 
count. Furthermore, the PRM treatment provides both 
fcrmionic and bosonic quasi-particle energies which are 
not directly accessible with other methods. However, the 
PRM approach of Ref. j2J was restricted to the metallic 
phase. 

Here, we extend our recent work |21| on the one-dimen- 
sional spinless Holstein model Q where the scope of this 
paper is twofold: Firstly, we demonstrate that the PRM 
approach of Ref. [2] is not restricted to the half-filled 
case, and that the phonon dispersion relation close to the 
critical value of the EP coupling reflects the type of the 
broken symmetry of the insulating phase. Secondly, we de- 
rive modified renormalization equations for the half-filled 
case that enable a dimerization of the system. Thus, we 
find an uniform description for the metallic as well as for 
the insulating phase of the spinless Holstein model QJ 
at half-filling. We determine the critical coupling g c for 
the metal insulator transition where a careful finite-size 
scaling is performed. Furthermore, a phonon softening is 
found for EP couplings close to the critical value g c which 
can be understood as a precursor effect of the phase tran- 
sition. 

This paper is organized as follows. In the next section 
we briefly describe the basic idea of our recently devel- 
oped PRM approach JH], and discuss metallic solutions 
of the renormalization equations for different fillings of the 
electronic band. In particular, we show that the phonon 
softening close to the critical value g c of the metal insula- 
tor transition reflects the type of the broken symmetry of 
the insulating phase. In sectional we extend the PRM ap- 
proach of Ref. |21j to derive an uniform description for the 
metallic as well as for the insulating phase of the spinless 
Holstein model Q at half-filling. Furthermore, in section 
the critical coupling g c for the metal insulator tran- 
sition is determined, and electronic as well as phononic 
quasi-particle energies arc presented. Finally, we summa- 
rize in section [3] 

2 Methodology and metallic solutions 

In the PRM approach |TB], the final effective Hamilto- 
nian TL = limA^o 7~Lx is obtained by a sequence of unitary 
transformations, 

n {X -AX) =e x ^H x e- x ^\ (2) 

by which transitions between eigenstates of the unper- 
turbed part TLq of the Hamiltonian caused by the inter- 
action TLx are eliminated in steps. The respective tran- 
sition energies are used as renormalization parameter A. 



The generator Xx.ax of the unitary transformation has to 
be adjusted in such a way so that TCrx-AX) does only con- 
tains excitations with energies smaller or equal (A — AX) . 
Interaction with energies larger than (A — AX) are used 
up to renormalize the parameters of the effective Hamil- 
tonian. Thus, difference equations for the A dependence of 
the parameters of the Hamiltonian can be derived which 
we call renormalization equations. 

In Ref. |2] we have made the following ansatz 

A = ^ £k,\c\c k + ^ Uq^bqbg + E\, (3) 
k q 

Wl,A = -7= 9k '1' X ( b l C l Ck +0 + b A+q C k) (4) 

k,q 

for the renormalized Hamiltonian TLx — TLq.x +7~ti \ after 
all excitations with energies larger than A have been elim- 
inated. In the next step, all excitations within the energy 
shell between (A — AX) and A have been removed where 

X X ,AX = ~7^J2 B MAZ1A (^4 C fc+9 - b q C l+q C k) 
k,q 

has been used as generator of the unitary transformation 
|J2J. The coefficients Bk, q ,x,AX have to be fixed in such a 
way so that only excitations with energies smaller than 
(A — AX) contribute to TCrx-AX)- 

Evaluating Eq. (J2J, operator terms that contain four 
fermionic and bosonic one-particle operators appear. How- 
ever, to restrict the renormalization scheme to the opera- 
tors of the renormalization ansatz according Eqs. J3J) and 
|1| we have to perform a factorization approximation, 

c ! c fc c L 9 c fc-9 ~ 4c fc (4_ 9 c fe _ 9 ) + (4c fc >4_ 9 cfc_ g 

-(4 C fc)(4-g C fc-9>> 

b\b q c{c k w b\b q {c\c k ) + {b\b q )c\c k - (b\b q ){ c {c k ) . 

Thus, the renormalization equations for £k,x, w ?,a, Ex, 
and gk,q,x, that are obtained by comparing the result of 
the transformation (J2J) with Eqs. © and J3J, depend on 
unknown expectation values (4 c fc) an d (b q b q ). These ex- 
pectation values are best defined with Hx since the renor- 
malization step is done from Tlx to Tin- AX) ■ Note that for 
simplicity in Ref. |21j the expectation values were evalu- 
ated with the unperturbed Hamiltonian 7io,A- In contrast, 
we now neglect the A dependence of the expectation values 
and perform the factorization approximation with the full 
Hamiltonian TL in order to take into account important in- 
teraction effects. (A general discussion of the factorization 
approximation in the PRM can be found in Ref. |22j.'l 

For the evaluation of the expectation values with the 
full Hamiltonian TL we have to apply the unitary transfor- 
mations also on operators to exploit (A) = {Ax)n x - The 
transformed operator Ax is obtained by the sequence (J2J) 
of unitary transformations, Aix-AX) = e. Xx ^ x Axe-~ Xx ^ x ■ 
We derive renormalization equations for the fermionic and 
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Fig. 1. (color online) Bosonic quasi-particle energies of a chain 
with 500 lattice sites for different EP couplings g obtained from 
the renormalization equations without a symmetry breaking 
term. Note the unphysical jumps at small q values which are 
due to the factorization approximation of higher order renor- 
malization processes. For details see the discussion in Ref. [21) . 

bosonic one-particle operators, c\ and b q , where the same 
approximations are used as for the Hamiltonian. In this 
way, equations for the needed expectation values are ob- 
tained. 

The resulting set of renormalization equations has to 
be solved self-consistently. The explicit (numerical) eval- 
uation starts from the cutoff A = A of the original model 
and proceeds down to A = 0. Note that the case A = 
with self-consistently determined expectation values pro- 
vides the effectively free model Tt — limA— o ~Ho,x, 

H = ^2e k clc k + ^2u q blb q + E, (5) 

k q 

we are interested in. Here, we have introduced the renor- 
malized dispersion relations e k = lirnx_>o£fc,A and uj q — 
lim>,^o w 9,a, and the energy shift E = lim.A-*o E\. 

Here, one might wonder why it is possible to map the 
Holstein model of spinless fermions onto an effectively free 
system as described above. Of course, before the actual 
calculations can be started a guess of the form of the final 
renormalizcd Hamiltonian TL is needed that can be moti- 
vated by the operator terms generated due to the renor- 
malization procedure. Furthermore, the obtained results 
are a very powerful test of the inital guess: Unphysical 
findings are clear signatures of an insufficient renormal- 
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Fig. 2. Critical value g c of the EP coupling as a function of 
the electronic band filling for a chain with 500 lattice sites. The 
value g c w 0.31t for the half-filled case agrees perfectly with 
the result of Ref. |21[ where a simpler scheme for the evaluation 
of the expectation values has been used. 



ization ansatz. As an example, it will turn out that non- 
physical negative phonon energies uj q are obtained if the 
electron-phonon coupling g exceeds a critical value. There- 
fore, we shall later extend ansatz (3), Q to study both 
the metallic and the insulating phase of model (JIJ. 

Furthermore, it is important to note that the employed 
factorization approximation is directly related to the re- 
normalization ansatz: Only operator structures that are 
contained in the renormalization ansatz TL\ can be taken 
into account without any approximation. All other op- 
erators appearing due to the renormalization procedure 
have to be traced back to those of Tt\. Here, as already 
mentioned, we use a factorization approximation for this 
purpose. Due to the complexity of the renormalization 
equations it is extremely difficult to estimate the effect 
of the mentioned factorization approximation on the re- 
sults. However, in this paper we present two complemen- 
tary renormalization schemes so that the comparison of 
the two approaches will provide some information about 
the effects of the employed factorization approximations 
(see the discussion in Sec.0J. 

In the following we concentrate on the so-called adi- 
abatic case u>o <C t, where we have chosen Wo/t = 0.1. 
In panel (a) of Fig. ^ t ne bosonic dispersion relation of 
an one-dimensional chain at half-filling is shown for dif- 
ferent EP couplings g. Due to the coupling between the 
phononic and electronic degrees of freedom, Cj q gains some 
dispersion. In particular, a phonon softening appears at 
the Brillouin-zone boundary if g is only slightly smaller 
than the critical value g c « 0.31£ of the EP coupling. 
However, if the EP coupling g exceeds the critical value g c 
we obtain non-physical negative phonon energies. There- 
fore, the solutions of the renormalization equations from 
Ref. |2] can only describe the metallic phase and break 
down at the phase transition. 
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As already mentioned above, in the present calcula- 
tion an improved scheme for the evaluation of the expec- 
tation values has been used. However, the same critical 
value g c w 0.31t for the metal-insulator transition was 
also found in Ref. [22 where all expectation values had 
been evaluated using the unperturbed part Ho,x of the A 
dependent Hamiltonian and not using the full Hamilto- 
nian TL. 

Note that the critical EP coupling g c as 0.31t ob- 
tained from the vanishing phonon phonon mode at the 
Brillouin-zone boundary is significantly larger than the 
DMRG value of g c w 0.28t [TDIH] . This difference can 
be understood as follows: Our renormalization scheme as 
presented in this Sec. starts from a description partic- 
ularly suitable for the metallic phase as represented by 
the final Hamiltonian H. As discussed above, a factoriza- 
tion approximation has been employed so that fluctua- 
tions are suppressed. However, these additional fluctua- 
tions would enhance the phonon softening and, therefore, 
tend to destabilize the metallic phase. Thus, the realistic 
critical value g c of the electron-phonon coupling is smaller 
than the result we have obtained by the PRM approach as 
presented above. Note that the situation will change for 
the modified renormalization scheme of Sec. |5] 

Even though only the case of half-filling has been dis- 
cussed in Ref. [21], the same PRM approach can also 
be applied to different fillings without any modifications. 
Panel (b) of Fig. shows bosonic dispersion relations for 
an one-dimensional chain at filling 1/3 for different EP 
couplings. Similar to the half-filled case, we obtain a pho- 
non softening for EP couplings g close to a critical value 
of g c ~ 0.5t, and negative phonon energies are observed 
for g > g c . However, in contrast to the half-filled case, 
the phonon softening appears at 2Uf = 2tt/3 and not 
at the Brillouin-zone boundary (where kp denotes the 
Fermi momentum). Consequently, as expected for a one- 
dimensional fermionic system, the broken symmetry of the 
insulating phase strongly depends on the filling of the elec- 
tronic band. 

Fig. El shows the critical values g c of the EP coupling 
as function of the electronic band filling. Here, two aspects 
of the results are important to be noticed: (i) The critical 
value g c has a minimum at half-filling which is connected 
with the highest stability of the insulating phase, (ii) The 
form of the curve in Fig. [21 reflects the particle-hole sym- 
metry of the model JTJ. 

Above we have argued that the phonon softening is 
a precursor effect of the metal-insulation transition in the 
Holstein model. However, in the anti-adiabatic limit of the 
model where t <C wo holds we do obtain the opposite be- 
havior: A phonon stiffening occurs in the metallic phase 
if we increase the electron-phonon coupling g. This find- 
ing corresponds to the anti-adiabatic limit of spin systems 
where also no phonon softening occurs for fast phonons 
|23| . Unfortunately, in this work the metal- insulator tran- 
sition itself can not be studied in the anti-adiabatic limit 
because the presented PRM approach has no stable solu- 
tions for t <C g. 



3 Uniform description of metallic and 
insulating phases at half-filling 

In this section we want to extend our PRM approach to 
find an uniform description for the metallic as well as for 
the insulating phase of the spinless Holstein model (TIJ. 
In the following, we concentrate on the case of half-filling 
because the broken symmetry of the insulating phase de- 
pends strongly on the filling of the electronic band as dis- 
cussed above. 

We have already mentioned that the PRM approach as 
discussed in Sec. [21 breaks down for strong EP couplings. 
In this case a long-range charge density wave order oc- 
curs, and electrons and ions are shifted from their sym- 
metric positions. To describe such a state with a broken 
symmetry in the framework of the PRM the ansatz for 
Tlx must contain suited order parameters |20| . In the case 
of the half-filled Holstein model, one has to take into ac- 
count that the unit cell of the system is doubled in the 
case of a dimerized insulating ground-state. Therefore, we 
consider the Hamiltonian in a reduced Brillouin zone and 
introduce appropriate symmetry breaking terms in our 
renormalization ansatz. Thus, the renormalized Hamilto- 
nian H\ = TLo,\ + Tii, a reads 

n ,x = (6) 

= ^2 £ n,k,\Ca,k C <*,k + E W 7>9>A & 7,g & 7,9 + E * 
fc>0,Q 9>0,7 

+ *%,x (4,^,* + h.c.) + VNAl (b{ Q + h.c.) 
k 

= 

= ~m E 9t: j{^ q )5{cl k c p , k+q ) + h.c.} (7) 

V fc,g>0 

after all excitations with energies larger than A have been 
integrated out. Here, due to the usage of the reduced Bril- 
louin zone, both the fermionic and bosonic one-particle 
operators as well as the model parameters have additional 
band indices, a, (3, 7 = 0,1. Furthermore, the definitions 
5 A = A - (A) and Q = n/a are used in Eq. 0. Note, 
however, that Eq. JHl is restricted to the one-dimensional 
case at half-filling but could be easily extended: For higher 
dimensions the term with the order parameter has to 
take into account all Q values of the Brillouin zone bound- 
ary. To describe the insulating phase away from half-filling 
one has to adjust the order parameters because the broken 
symmetry type strongly depends on the filling as discussed 
above. 

To derive the renormalization equations we consider 
the renormalization step from A to (A — AX). At first we 
perform a rotation in the fermionic subspace and a trans- 
lation to new ionic equilibrium positions so that Tlo,x is 
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diagonalized, 

Hq.X = ^2^2£a,k^ C i,k,X C a,k,\ (8) 

fe>0 a 

q>0 7 

Next we rewrite TL\^\ in terms of the new fcrmionic and 
bosonic creation and annihilation operators, C^ k > and 

A which we have introduced to diagonalize 7io,A- Fi- 
nally, we have to evaluate to derive the renormalization 
equations for the parameters of H x . Here the ansatz 

k,q a,/3,7 

x{sB\ <q 5(Cl x C 0tk+q , x )-h.c.} 

is used for the generator of the unitary transformation © . 
The coefficients A k '^' x Ax have to be fixed in such a way so 
that only excitations with energies smaller than (A — AX) 
contribute to 7~ii.(x-AX)- The renormalization equations 
for the parameters e k .x, A k A , ui 7 . q _\, A b x , and g kq x are 
obtained by comparison with ©, after the creation 
and annihilation operators C^ k A , B^ A have been trans- 
formed back to the original operators c£ fc , b^} q . The actual 
calculations are done in close analogy to Ref. [21] ■ 

As already discussed in Sec. [21 i n order to evaluate 
Eq. (J2J, an additional factorization approximation must 
be employed where only operators of the structure of those 
of © and (0) are kept. At this point it is important to no- 
tice that this factorization approximation is performed in 
the framework of the reduced Brillouin zone. Therefore, in 
comparison to the PRM scheme of Sec. [21 additional terms 
occur for non-zero order parameters A k x and A\ . Due to 
the factorization approximation, the final renormalization 
equations still depend on unknown expectation values. As 
already discussed in Sec. [21 these expectation values were 
evaluated with the full Hamiltonian TL in order to take 
into account important interaction effects. Therefore, we 
again must apply the sequence 121 of unitary transfor- 
mations to operators, A(x-AX) = e Xx - AX Axe~ XxAX to 
exploit (.4) = (Ax)hx- As in Sec. [21 this procedure is 
performed for the fermionic and bosonic one-particle op- 
erators, c' a k and &t , where the same approximations are 
used as for the Hamiltonian. Thus, we easily obtain equa- 
tions for the needed expectation values. The resulting set 
of renormalization equations is solved numerically where 
the equations for the expectation values are taken into 
account due to a self-consistency loop. 

The case A = with self-consistently determined ex- 
pectation values provides again an effectively free model 



H = lim^-jo Hx = liniA^o ^o,a which reads 

h = E ^,k4,, k ca,k + J2 A i (4,fe c i,fe + h - c -) ( 9 ) 

k>0,a k>0 

+ "iAaK* + (b\ Q + & 1)Q ) - E 

9>0,7 

where e a , k = lim A ^o e a ,k, x, A c k = lim A ^o^A> ^i,q = 
lim;\_>o UJ 1<q< \, and A b = limA^o A\. Note that all excita- 
tions from TL\ t x were used up to renormalize the param- 
eters of 7io- The expectation values are also calculated 
in the limit A — > and can be easily determined from 
(A) H = {Ax)h x = ((lim A ^o because H is a free 

model. 

Before we present results in Sec. 01 we want to com- 
pare the two different renormalization schemes discussed 
above in more detail. As already mentioned, the approach 
of Sec.[2]is based on a renormalization ansatz in particular 
suitable for the metallic phase of the system. In contrast, 
the renormalization scheme as discussed in the present 
section starts from a dimerized system so that this ap- 
proach provides a description of the system particularly 
adapted to the insulating phase. Thus, both renormaliza- 
tion schemes are complementary approaches, and opposite 
fluctuations are suppressed due to the employed factoriza- 
tion approximations: The suppressed fluctuations of the 
renormalization scheme of Sec. [21 destabilize the metallic 
phase. In contrast, it will be shown that the neglected fluc- 
tuations of the approach of Sec. 01 favor the same phase. 

On the other hand, the two renormalization schemes 
are not only complementary in their starting points but 
also closely related: The renormalization scheme as dis- 
cussed in the present section exactly matches the approach 
of Sec. [2]if the gap parameters Ai x and A x vanish for all 
A values. Thus, finite excitation gaps are only obtained 
from the renormalization scheme of this section if the in- 
sulating (i.e. gapped) phase has a lower free energy than 
the metallic solution as discussed in Sec. [21 Note that in 
actual calculations the PRM approach of this section al- 
ways leads to a free energy smaller than the one obtained 
from the scheme discussed in Sec. [21 



4 Results 

In the following, we first show that the PRM treatment 
can be used to investigate the Peierls transition of the one- 
dimensional spinless Holstein model at half-filling. In 
particular, our analytical approach provides a theoretical 
description for the metallic as well as for the insulating 
phase, and we compare our results with DMRG calcula- 
tions [TTllB] . 

In order to determine g c we perform a careful finite-size 
scaling as shown for some g values in the inset of Fig. 
where a linear regression is applied to extrapolate our re- 
sults to infinite system size. Note, however, that the finite 
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Fig. 3. Gap in the electronic excitation spectrum for the in- 
finite chain. The solid line is a Kosterlitz-Thouless, fit of the 

15.240 



form A/t 



V( 9 /t)2-(0.195)' 



exp 



1.400 



V( 9 /t)2-(0.195) 2 



The in- 



set shows the finite-size scaling for g values of 0.26i (circles), 
0.2664 (squares), and 0.274 (diamonds). 



size scaling of our PRM approach is affected by two differ- 
ent effects: Long-range fluctuations are suppressed by the 
finite cluster size as well as by the used factorization ap- 
proximation so that a rather unusual dependence on the 
system size is found. In Fig. [3] the electronic excitation 
gap A for infinite system size, as found from the opening 
of a gap in the quasi-particle energy ik (see text below), 
is plotted as function of the EP coupling g. A closer in- 
spection of the data shows that an insulating phase with 
a finite excitation gap is obtained for g values larger than 
the critical EP coupling g c sw 0.244. A comparison with 
the critical value g c « 0.284 obtained from DMRG cal- 
culations shows that the critical values from the 
PRM approach of Sec. El might be somewhat too small. 

In particular, the critical coupling g c sa 0.244 obtained 
from the opening of the gap in Sk is significantly smaller 
than the g c value of ~ 0.314 which was found from the van- 
ishing of the phonon mode at the Brillouin zone boundary 
in the metallic solution of Sect. [51 Instead, one would ex- 
pect that both the gap in ik and the vanishing of uj q should 
occur at the same g c value. This inconsistency can be un- 
derstood from the different nature of the factorization ap- 
proximations employed in the two PRM approaches pre- 
sented here: As discussed above, the due to the factoriza- 
tion approximation suppressed fluctuations tend to desta- 
bilize the suggested phase of the renormalization ansatz. 
Therefore, the neglected fluctuations of the two presented 
PRM approaches have opposite characters: The stability 
of the metallic phase is overestimated by the scheme of 
Sec. |2 whereas the approach of Sec. |3 favors the insulat- 
ing phase, and the different values of the critical coupling 
exactly reflect the different nature of the two PRM ap- 
proaches. Thus, both approaches (and both ways to deter- 
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Fig. 4. (color online) Fermionic (upper panel) and bosonic 
(lower panel) quasi-particle energies of a chain with 500 lattice 
sites for different EP couplings g. 



mine g c ) would be consistent with each other if additional 
fluctuation could be taken into account, and one would 
expect a common result for g c between 0.244 and 0.314. 
This would be in good agreement with the DMRG value 
of g c « 0.284 HHEj. 

Another way to determine g c might be given by the as- 
sumption that a marginal relevant interaction in a renor- 
malization group treatment causes the phase transition in 
the Holstein model. Corresponding to Ref. (21), the gap A 
should in this case grow for g > g c as exp[— a/g— g c ] , where 
a is some numerical constant. Such a function fits our data 
very well, and a critical value of g c = 0.1664 would result 
in this way which seems to be questionable small. On the 
other hand, we could also assume that the metal-insulator 
transition might be of Kosterlitz-Thouless type [2SJ as 
found for the anti-adiabatic limit of spin-Peilerls chains 
|26| . As one can see in Fig.0 the Kosterlitz-Thouless gap 
formula |27ll2l>| .^(x (g 2 - 



-9\ 



') - 5 exp(-a/^ 5 2 - ff 2) ;a i so 
fits our data and leads to a critical value of g c = 0.1954. 
However, such a small value contradicts other finding of 
our calculations: As discussed above, the phonon softening 
at the Brillouin-zone boundary is a clear signature of the 
occuring phase transition, and the smallest phonon ener- 
gies should be obtained for g c . This test provides clear 
evidence for a critical value of g c « 0.244. 

The quasi-particle energies are also directly accessi- 
ble: After the renormalization equations were solved self- 
consistently the electronic and phononic quasi-particle en- 
ergies of the system, ik and u> q , respectively, are given by 
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the limit A — » of the parameters £a,k,x an d w x9,> °^ ^ e 
unperturbed part 7io,A of the Hamiltonian in its diagonal- 
ized form (JSJ). 

In Fig.rjjthe renormalized one-particle energies (which 
have to be interpreted as the quasi-particle energies of the 
full system) are shown for different values of the EP cou- 
pling g. As one can see from the upper panel, the electronic 
one-particle energies depend only slightly on g as long as g 
is smaller than the critical value g c ss 0.24i. If the EP cou- 
pling g is further increased a gap A opens at the Fermi en- 
ergy so that the system becomes an insulator. Remember 
that the gap A has been used as order parameter to de- 
termine the critical EP coupling g c of the metal-insulator 
transition (see Fig. O ■ The lower panel of Fig. [J] shows the 
results for the phononic one-particle energies uj q . One can 
see that u> q gains dispersion due to the coupling g between 
the electronic and phononic degrees of freedom. In partic- 
ular, the phonon mode at momentum 2kp = 7r, i.e. at the 
Brillouin-zone boundary, becomes soft if the EP coupling 
is increased up to g c « 0.24t. However, in contrast to the 
metallic solution of Sect. at 2kp always remains pos- 
itive though it is very small. Note that for g values larger 
than g c Cb q increases again. This phonon softening at the 
phase transition has to be interpreted as a lattice instabil- 
ity which leads to the formation of the insulating Peierls 
state for g > g c . The phase transition is associated with a 
shift of the ionic equilibrium positions. A lattice stiffening 
occurs if g is further increased to values much larger than 
the critical value g c ~ 0.24t. 

5 Summary 

In this paper, the recently developed PRM approach has 
been applied to the one-dimensional Holstein model of 
spinlcss fcrmions interacting with dispersion-less phonons. 
In extension to our earlier work |21|. here, we have im- 
proved the scheme for the evaluation of expectation val- 
ues, and have discussed the metal-insulator transition of 
the model for different fillings of the electronic band. Fur- 
thermore, for the half-filled Holstein model we have de- 
rived an uniform description for the metallic as well as for 
the insulating phase. 

We have shown that the renormalized phonon energies 
gain a momentum dependence due to the EP coupling. In 
particular, for half-filling a phonon softening at 2kp ap- 
pears if the EP coupling is close to the critical value of the 
metal-insulator transition. Therefore, the broken symme- 
try of the insulating phase strongly depends on the filling 
of the electronic band. The critical value of the metal- 
insulator transition also depends on the band filling, where 
the insulating phase has at half-filling the highest stability. 

The PRM approach of Ref. |2T] breaks down if the EP 
coupling exceeds the critical value of the metal-insulator 
transition. Here, we have extended our PRM approach to 
enable symmetry broken insulating states. Thus, we de- 
rived an uniform description for the metal as well as for 
the insulating phase (which is restricted to the half-filled 
case of the Holstein model). Here, we have used this ex- 
tended analytical framework to study the phase transition 



in more detail. We have determined the critical value of 
the EP coupling for the metal-insulator transition. Fur- 
thermore, we have shown the opening of the excitation 
gap in the electronic quasi-particle energies, and the lat- 
tice stiffening for EP couplings larger than the critical 
value. 
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